Devil's staircase of incompressible electron states in a nanotube 
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It is shown that a periodic potential applied to a nanotube can lock electrons into incompressible 
states. Depending on whether electrons are weakly or tightly bound to the potential, excitation gaps 
open up either due to the Bragg diffraction enhanced by the Tomonaga - Luttinger correlations, or 
via pinning of the Wigner crystal. Incompressible states can be detected in a Thouless pump setup, 
in which a slowly moving periodic potential induces quantized current, with a possibility to pump 
on average a fraction of an electron per cycle as a result of interactions. 
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A carbon nanotube (NT) is a strongly interacting elec- 
tron system known to be a host of many-body effects [1] , 
including Luttinger liquid behavior [2] . At half- filling the 
NT is effectively dilute, and a crossover from the liquid 
to the Id Wigner crystal [3] is expected [4] . 

Here we suggest that an external periodic potential can 
be a probe of both crystallization and Luttinger correla- 
tions. We show that incompressible electron states arise 
when the electron number density p (relative to half- 
filling) is commensurate with the potential period Aoxt ■ 



P = "T-tot/Acxt , "T-tot = 4m . 



(1) 



In Eq. (1), m is the number of the NT electrons of each 
of the four polarizations [5] per period. To calculate exci- 
tation gaps we generalize the Pokrovsky-Talapov theory 
[6] for the case of the four coupled fermion modes. 

In the absence of interactions, Bragg diffraction on the 
potential opens minigaps at integer density (1), |m| = 
1,2, ... [7]. Interactions profoundly change the spectrum, 
yielding a devil's staircase of incompressible states at ra- 
tional m = p/q. In such a state, the NT electron system 
is locked by the potential into a gAoxt^periodic structure. 
If detected, e.g. in a Thouless pump setup [7], corre- 
sponding minigaps would provide a direct probe of inter- 
actions, with a possibility to map the devil's staircase by 
pumping at fractions of the base frequency. 

One-dimensional interacting electrons are convention- 
ally described by the Tomonaga - Luttinger liquid [2]. 
This hydrodynamic approach is valid in a small momen- 
tum shell near the Fermi points, with excitations ex- 
tended over the whole system. Adequate description of 
crystallization and commensurability requires including 
the curvature of the electronic dispersion that becomes 
important at low density. The curvature can yield crys- 
tallization or commensuration by coupling charge and 
spin modes and by introducing a length scale into an oth- 
erwise scale-invariant Gaussian theory. In this work we 
treat both electron interactions and the curvature of the 
dispersion non-perturbatively by making use of the rel- 
ativistic Dirac spectrum of a half-filled nanotube. Cur- 
vature is controlled by the Dirac gap and is bosonized 
exactly by virtue of the massive Thirring - sine-Gordon 



duality. [8, 9] This enables us to study incompressible 
states both in the limit of a narrow-gap Luttinger liquid 
and in that of the locked Wigner crystal. 

Our course of action is to introduce the bosonized de- 
scription for the NT electrons, develop the phase soliton 
method and find excitation gaps from the renormalized 
sine-Gordon action, draw the phase diagram in the semi- 
classical limit, and comment on the experimental means 
to detect incompressible states. 

The model. — Nanotube electrons in the forward scat- 
tering approximation [10] are described by the four fla- 
vors of Dirac fermions V'a = (V'q V'a) , a = 1 — 4, whose 
interaction is written in terms of the smooth envelope 
P(^) — Sl=i '/'^ (a^)V'ct(a;) of the total charge density. 
The second-quantized Hamiltonian Ti = Wo+'^bs + Wcxt, 
where Tig is the massless Dirac Hamiltonian 



TLq = -ihv I ^ ^^asdxipadx + ^ X! P-kV{k)pk (2) 



with the Goulomb interaction V{k) — In [l + (/ca)"^] 
for a tube of radius a placed on a substrate with the di- 
electric constant e. The curvature of the electron disper- 
sion controlled by the gap 2Ao at half-filling introduces 
backscattering at each NT Dirac point: 



(3) 



We emphasize that the backscattering Aq is not the usual 
interaction-induced V{2kF) term (which is undetectably 
small in metallic NTs), but rather is present at the single- 
particle level [5, 11]: Depending on the tube chirality, 
the bare gap Aq can be in the range Aq <10meV to 
'--^O.S eV; it can also be controlled by the parallel magnetic 
field [12]. Adding the periodic potential U{x) and the 
chemical potential p [p = at half-filling] results in 



U{x) 



Hc^t = J dxp{x){U{x)-p} , (4) 
= A cos ko^tx , fccxt = Stt/ Aoxt • (5) 



Qualitatively, our findings will be valid for any realistic 
potential which justifies the simplest choice (5); typically. 
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Aext 0-1 ^ 1 MHi- The Hamiltonian Ti. is U(4) invariant 
with respect to rotations in the fermion flavor space. 

Bosonization of the nanotube electrons Va = 
--^== e'®" is exact even in the presence of (3). It maps the 
problem of the four interacting Dirac fermion modes onto 
the sine-Gordon model of the four coupled bose fields ©„. 
We rotate [4, 7] to the basis of the charge mode 0° and 
three neutral modes O"', in which case the charge density 
p{x) = ^ dxO^, and the Gaussian action (bosonized Ho) 
is diagonal [Ji = v = 1]: 

jCo = ^Jdx [{d,9'r - K {djy + ELi {9,0"?) ■ 

(6) 

The slow momentum dependence of the charge stiffness 
Kk = 1 + AvV{k), = nhv, is irrelevant, and we 
take it as constant K = K^^^iji^^, Zch being the size of 
the charged- mode soliton (described below). Assuming 
^ch ~ Is, with the screening length Zg ~ 1 /xm, and using 

/hv ~ 2.7, one estimates K AQ for for the stand- 
alone tube; ii" ~ 10 if the tube is placed on a substrate 
with dielectric constant £ ~ 10. The logarithmic behav- 
ior of V{k) underlies the fact that the Coulomb interac- 
tion V{x?) oc l/\x\ is essentially local in one dimension, 
\ J dxdx'p{x)V{x - x')p{x') / dxK{dJ°f. 

The nonlinear part of the sine-Gordon Lagrangian 
[coming from Wbs] reads [4, 7, 13] 

£bs = - j dxgoT (e° + 2pk^^tx -2Asmk^^tX, ,(7) 

JP(6'°, =cos6'"-na=iCos6'" + sin0O-nLiSin6'" ,(8) 

where go = AAo/irDa^ , and D c± hv/ais the Id band- 
width. In Eq. (7) we included the coupling (4) to ex- 
ternal fields by shifting the charge mode 9^ ^ 9^ — 
^fK-\U- p) dx', with p = p/{Keo), A = A/{Keo), 
eo = Tifcextf . 

In what follows, it is useful to first describe elementary 
excitations of the stand-alone tube, U = Q. As shown by 
Levitov and Tsvelik [4] , in the bosonized language adding 
one electron corresponds to a com,posite soliton of both 
the charge and the fiavor modes. In such a composite 
object, the charge mode changes by 7r/2 (adding unit 
charge) over the length Zchj whereas the neutral sector 
adds a particular SU(4) fiavor to the electron by means 
of the solitons of 9°" "switching" by ±7r/2 on a shorter 
scale l[\ ^ K~^/'^lc\i, right in the middle of the charge 
soliton. The composite soliton is a unit charge configu- 
ration of the minimal energy, obtained by optimizing the 
action Co + C^s m the limit of large Coulomb repulsion 
K ^ 1. Technically, such an optimization results in the 
soft neutral modes adjusting to create the effective 
potential T{9^) = io.m{0a} J^{9°- , 9°) ~ cos46l° for the 
stiff charge mode 9°. 

In the absence of external potential, U = 0, the system 
Co + Chs describes the Wigner crystal - Luttinger liquid 
crossover. In particular, raising the chemical potential 



from half- filling [band insulator, or the "Dirac vacuum"] 
to just above the charge gap produces the train of weakly 
overlapping composite solitons described above [pla <^ 1, 
"Wigner crystal"], in which overlapping charge solitons 
maintain a quasi-long-range order. Further increase of p 
leads to strongly overlapping solitons rendering the non- 
linear term (7) irrelevant [pZfl ;S> 1, Luttinger liquid]. We 
emphasize that the crystal - liquid crossover occurs due 
to the finite soliton size oc Aq ^ that scales inversely 
with the curvature of the dispersion. Curvature is also 
responsible for binding flavor to charge through the term 
(7). The U = system is compressible [4]. 

Periodic potential locks electrons into incompressible 
states. Technically, the term (7) becomes relevant when- 
ever the density 2p, = m [Eq. (1)] is integer [7]. When 
m is a simple fraction, 2/i = p/q, to identify incom- 
pressible states one utilizes the phase soliton method 
[6]. We generalize it for the case of the four nanotube 
modes by expanding in powers of the coupling go 
6»J = 9^ + go9^^'^^ + ... + gp^^^^l j = 0, a, and flnding 
the effective Lagrang ian Cm \9^ , 9^\ for the phase modes 
0^ and 9°' , a = 1—3, to the lowest order in go- The 
phase modes are constant in the commensurate 

phase, whereas an excitation is a composite phase soli- 
ton, in which the phase fields ^"(x), 6"'(x) describe a slow 
deformation of the regular commensurate configuration. 
Excitation gap is given by the energy of the composite 
soliton, renormalized by quantum fluctuations. 

We note that it is the curvature cx Aq that yields in- 
compressible states. Same is true for the non-interacting 
electrons: When Aq = 0, the external potential is gauged 
away from the Hamiltonian —ikva^dx + U{x). 

The composite phase soliton is a result of optimization 
of the corresponding effective action Cm[9'^,0°'] (exam- 
ples of which are given below). This problem will be 
similar to the U = case described above: When in- 
teraction is strong. A" ^ 1, the neutral phase modes 0" 
adjust [on the scale la] to create an effective potential 
for the charged phase mode 9^. The crucial difference 
from the former case is that this optimization will quali- 
tatively depend on whether the system is in the regime of 
the Luttinger liquid or in that of the Wigner crystal. In 
this sense, periodic potential naturally distinguishes be- 
tween the opposite sides of the crystal - liquid crossover, 
by bringing about the additional length scale, its period 
Aext- Technically, the saddle point of Cm will depend on 
whether the neutral modes are adiabatic or fast on the 
length scale on which the effective potential (7) changes 
appreciably. Below we consider both cases separately. 

Bragg diffraction in a four-flavor Luttinger liquid. — 
In the adiabatic limit la » Aext of extended flavor exci- 
tations, exchange is important: The system (correlated 
over many Aext) "knows" that it is comprised of particles 
of the four different flavors. This limit is connected to 
the non-interacting case, where particles repel only due 
to the Pauli principle. For integer density m, averag- 
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ing (7) over the period Aext reduces the problem to the 
U = case with = + ^hs, go goJ-mC^A), yield- 
ing the » 1 minigaps A„ ~ K^/^D^/^\A^^^\^/^ that 
are enhanced by the quantum fluctuations compared to 
the noninteracting values = 2AoJto(2A) [7]. 

The fractional m = | case is considered in de- 
tail in Ref. 13. A somewhat lengthy calculation 
yields the effective action of the form £1/2 = £0 + 
v{A) J2a COS 29'' cos 20" - u{A) J2a>b COS 29" cos 29''. We 
find the charge excitation gap Ax/2 oc D 
(that vanishes in the noninteracting limit /sT — > 1 in 
accord with the Bloch theory), whereas flavor excita- 
tions arc governed by the SU(4)~0(6) Gross - Neveu 
Lagrangian [14] that can be written in terms of six 
Majorana fermions x ~ e'^, £gn = iXjlf^df^Xj - 
'^{■^){XjXj){X3'Xj')- Although the charge gap A1/2 van- 
ishes for certain values of the potential amplitude (zeroes 
of v{A) given in Ref. 13), the excitation gap never closes 
due to flavor, since the Gross - Neveu coupling u{A) ^ 0. 

Locked Wigner crystal. — Below wc focus on the op- 
posite limit la -C Aext, in which the charge mode domi- 
nates, whereas exchange (overlap of the flavor solitons) 
is relatively unimportant. Wc will show that the effec- 
tive description in this regime is that of a single mode of 
fermions with the density (1), locked by the external po- 
tential. Physically this happens since when the electron 
wavefunctions (represented by solitons) are localized on 
the scale Aext, Coulomb repulsion wins over exchange. 
Naturally, in this case the appropriate saddle point of 
Co + Cbs is equivalent to the semi-classical description, 
with the WKB condition Zfl <C Aext- Technically, the 
neutral mode "switching" produces the effective charge 
mode potential with a quarter-period, similarly to the 
U = case above: gfo-^(^'^+"^fcexta^ — 2AsinA;ext2;, ^") ~ 
g cos(46''' + mtotkcxtx — SA sin fcoxta;) . This potential now 
depends on the total density mtot, Eq- (!)• 

To obtain the renormalized coupling g we utilize adia- 
baticity of U{x), la <C Aext- The renormalization group 
produces the effective coiipling g on the scales a < I < Ifi 
that appear "microscopic" for the external potential. 
Thus the neutral modes can be integrated out, yielding 
both Ifi and g independent of the shape of the potential 
U{x). The problem becomes similar to the U = case. 
The coupling g follows from the scaling g ~ go{lfi/o,)~^^^ 
and self-consistency gQa) ^ 1/Zfl'. [4, 7, 13] 

Re-fermionization. — Let us now complete our single- 
mode description by mapping the system in the limit 
Zfl <C Aext onto the problem of spinless Dirac fermions 
in the external potential (5). Wc introduce the displace- 
ment field e = 29° for the total density p = ^d^Q. To 
preserve commutation relations, we rescale the canonical 
momentum He = jHeo. Changing variables in the La- 
grangian 2^ {K{dx9°)^ + gcos{4:9° + ...)}, we obtain the 



effective Lagrangian for the charge mode 9, 

£,ff[9] = ^ / dx\-^{dtef ^ {d.&f 

-fif cos(29 + mtotke^tx - 2A' sin fcext^;) > 

with the rescaled parameters v' = 4u, K' = K/16, and 
A' = A/ K'Tikcxiv' = AA. The velocity quadrupling 
corresponds to counting incoming fermions regardless of 
their (four) flavors. With the same accuracy that al- 
lowed us to neglect the flavor sector, the rescaled value 
K' 1 + {K -1)/16= 1 + v'V{q) is by definition the 
charge stiffness for the spinless Dirac fermions of velocity 
v' , with the density of states v' = \/tt%v' . The external 
fields U and ji are not rescaled. Introducing the Dirac 
spinors = ,} , e^^ relative to the new cutoff a' /g, 
we obtain the effective Hamiltonian 

Weff [*] = fdx^+ {-ihv'aad^ + AVi + U{x) 

1 ("^^) 

k 

for the fictitious spinless Dirac fermions of unit charge, 

with the number density p{x) = coupled to the 

external fields in the natural way, Eq. (4). The fermions 
^' represent the original NT electrons traced over the 
flavor states. The effective gap A' in Eq. (10) corresponds 
to the renormalized coupling g entering the Lagrangian 
(9) hygc^^, yielding A' ^ ^ DV5A*/^ 

The (charge) excitation gaps are now estimated via 
the phase soliton approach [6] for the single-mode system 
(9) and (10), yielding the devil's staircase. For integer 
mtot = ±1, ±2, in the classical limit — > 00, we find 

A^:.K'/Vm.^.W\'/'D'/'A'J' . (11) 

The effective density quadrupling m mtot in (H) is 
due to strong interactions: Electrons of all flavors avoid 
each other, their ground state being a Slater determinant 
of the quadrupled size. 

Phase diagram. — The Hamiltonian (10) in the limit 
K ^ 1, coarse-grained beyond la, describes the clas- 
sical chain of rig electrons and Uh holes per period, 
'T^tot = rie — Uh- Their positions Xi in the minima and 
yj in the maxima of U (x) , as well as the characteristic 
Coulomb blockade hexagons labeled by {ue, rih) in Fig. 1, 
are obtained by minimizing Eci{A)/N—iJ,{ne—nh), where 

E,, = J2[^' + U{x.)] + ^[A'-[/(y,)] 

+ J2 ^(^* - ^«') + E ^(y^ - yr) - E ^(^^ - y^) ^^^^ 

i>i' j>j' i,j 

for a finite tube of the length L = NXcxt- In Fig. 1, the 
hexagon size is controlled by the NT capacitance per pe- 
riod, Co — 2in(i°]x t) ' n-e,nh ^ 1, the Thomas - 
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FIG. 1: Left: Phase diagram in the {A, /i) plane according 
to the model (12). In the limit ne,nh 3> 1 the period in A 
becomes e^/Ci. Right: Mapping the devil's staircase in the 
Thouless pump setup. In addition to the integer-m plateaus, 
the ones with fractional m appear due to electron interactions. 



Fermi approximation p{x) w Ptf{x) oc U{x) estimates 
the interaction energy inside each half-period ("quan- 
tum dot") as 4 fdxdx' ptfWptf(^') = ^here the 

"dot capacitance" Ci ~ 27r in(A''*t/a') • ^^^^ Hmit, all 
the hexagon domains of the phase diagram become iden- 
tical with the width e^/Ci. Remarkably, this period 
asymptotically follows from Eq. (11) with Jl' « ^ • ^sjo^ 
[K = i^i/A^^t ~ K — 1], using the Bessel function ze- 
roes 2i'^"^ 



+ -h Trn, n = min{ne, Uh}. The 
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mtot ' ~ 4 

borders between the domains correspond to the incom- 
pressible states with fractional n,, and Uh (not shown). 

Quantized current. — The devil's staircase can be 
mapped in the Thouless pump setup [7] . With the Fermi 
level in the gap, a slowly moving wave U{x — st) with the 
frequency / = s/Aext will generate the quantized current 
j = TTitotef- The moving potential can be created by 
gating, optical methods, or acoustic field. Estimated gap 
values A„ in the meV range ensure adiabaticity leading 
to current quantization and possible metrological appli- 
cations. Novel incompressible states with fractional m 
would correspond to the current quantized in the frac- 
tions of Aef, as illustrated in Fig. 1. Current changes 
sign at half-filling due to the Dirac symmetry. 

Conclusions. — We demonstrated that coupling to a 
periodic potential results in the devil's staircase of in- 
compressible electron states. Excitation gaps are found 
in the limit of the narrow-gap Luttinger liquid and in 
that of the Wigner crystal, by selecting the saddle point 
of the nonlinear action of the four bosonic modes. When 
the Coulomb interaction dominates, the system behaves 
as a single fermion mode with renormalized mass and ve- 
locity. Control over the NT gap, the screening length, 
and the parameters of the potential makes this setup a 



unique probe of the Luttinger liquid - Wigner crystal 
crossover and of commensuration effects in Id. Novel ef- 
fect of adiabatic pumping at fractions of base frequency 
is linked to the interaction-induced incompressible states. 
The effective single-mode description could also ratio- 
nalize recent manifestations of Wigner crystallization in 
transport, such as the e^//i steps in conductance [15]. 

This work was supported by NSF MRSEC grant DMR 
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